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ABSTRACT 

A study has been made on the possibility of identifying faulty 
components of a network by making tests only on external points of 
that network. The conclusion is made that this is possible, and that 
a practical method of doing this can be developed from the computer 
program presented as a result of this investigation. 

The method used is to select the voltage transfer functions of 
a network as the quantities on which the tests will be made. The 
poles and zeros of these functions are used to select a set of test 
frequencies. From the measurements made at these frequencies, 
a set of signatures is available which allow the faulty components 
to be identified. 
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L INTRODUCTION 


One of the urgent needs facing today^s Navy is a solution to the 
problem of electric and electronic equipment failures. It is not 
likely that such highly reliable components will be produced that the 
problem can be eliminated. Therefore, another approach is necessary. 

In the event of an equipment failure in an emergency situation, 
many lives may depend on how fast the repairs are made. Often 
the only method available to the technician is to disassemble the 

r 

equipment and make tests at various internal points. When the 
technician does not have a good familiarity with the equipment, 
this can lead to a down time of up to several hours, which in most 
cases is intolerable. Therefore, a method of maintenance is needed 
which would allow even an inexperienced technician to have the re^ 
pairs completed in the shortest possible time. Also, it would be 
desirable to be able to apply this same procedure to a routine 
program of preventive maintenance to eliminate these feilures before 
they can occur. 

A logical approach to this problem would be to devise a scheme 
that would use only external test points of the equipment -{i, e., the 
input and/or output terminals). Such an approach has been used in 
this investigation. Some previous studies in this area have been 
made [1], [ 2]. 
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This investigation has been carried out in the form of a feasi¬ 
bility study rather than to develop a practical method. The approach 
used was to select the voltage transfer functions of a network as 
the quantities on which tests are to be made. The considerations 
leading to this decision are given in Chapter II. 

Chapters III and IV cover the method of selecting test frequencies 
and a method of assigning signatures to the results obtained from 
these tests. Chapter V and the COMPUTER OUTPUT section 
cover the results of this investigation. 
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U. PRELIMINARY CONSIDERATIONS 


A. GENERATION OF NETWORK FUNCTIONS 

One of the initial steps undertaken in this study was the selection 
of a camputer program that generates netwark functions. The 
program chosen was giVen by Chang [3]. This program,^ written 
in the FORTRAN IV language, is based on topological formulas 
and uses the "dichotomy method" of a tree finding to generate the 
network functions. The program was modified to have it gene rate 
only the voltage transfer functions. 

The program to be presented as a result of this investigation 
did not substantially alter the tree-finding sections of Chang’s 
program. Therefore, a discussion of the theory and operation 
of the program taken from [3] will be given only when it is pertinent 
to the contents of this paper. 

B. PROGRAM LIMITATIONS 

Some limitations of the program,^ which will be presented, 
are now given. The program limits the*, networks that will be ad-^ 
mitted to those with graphs having the following properties: 

(a) non-oriented, (b) linear, (c) a maximum ©f 15 nodes, and: (d) 
a maximum of 60 edges [3]. 

The network functions are not generated in symbolic form 
(i. e.., the program requires that the inpuf data, the network elements,. 
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have numerical values). This requires that the tree-finding sub¬ 
routine be repeated each time that an element is varied, during 
the process of obtaining the data necessary for fault identification. 
The result of this is a substantial increase in the computer run 
time. 

These limitations do not impair the validity of this investiga¬ 
tion, since it is in the nature of a feasibility study. 

Additional comments on this subject will be made in a following 
section. 

C. CATEGORY OF NETWORKS CONSIDERED 

The limitations given in the preceding section place a restric¬ 
tion on the class of networks that will be allowed. The networks 
considered will be passive, linear networks (i. e. , containing only 
linear resistive, inductive, and capacitive elements). However, 
in this class of networks, the purely resistive networks will be 
prohibited, because of the frequency-dependent nature that is in¬ 
herent in the proposed test procedure. 

The assumption is also made that the networks under considera¬ 
tion have a defined set of input and output terminals. Thus, the 
permitted class of networks is: the set of passive linear, two-port 
networks. 

An illustration of a general two-port network is shown in 
Figure 1. The input and output terminals are represented by (1, 3) 
and (2, 4) respectively. 
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Figure 1 


The properties of two-port networks are well covered in the 
literature [3] - [6], The only two-port properties that are neces¬ 
sary for this investigation are the forward and reverse voltage- 
transfer functions. These are defined by equations (1) and (2) 
respectively. 

r - 

^11 - V 2 

D. VOLTAGE TRANSFER FUNCTIONS 

In the preceding section, it was stated that the only two-port 
properties needed were the voltage transfer functions. A nixmber 
of factors were considered before this decision was reached. 

First, the selection of the voltage transfer functions leads to 
a more simplified and convenient test procedure than would be ob¬ 
tained if some of the other network functions were selected. The 
test procedure proposed under this method will require only a 
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signal generator and a voltmeter. Obviously, using a combination 
of the other network functions, such as the z-parameters or cur- 
rent gains, would require more test equipment. Also, in the event 
of an emergency equipment failure, the test instruments most 
likely to be available would be a VTVM and a signal generator. 

This would be the only test equipment necessary in a program of 
preventive maintenance, which is the desired goal to be obtained 
from an extension of the test procedure presented here. The logi¬ 
cal conclusion reached from these points is to use the voltage trans¬ 
fer functions. 

In the preceding discussion, the suggestion was made to use 
the voltage transfer functions in the test procedure. The choice 
must then be made to use one or both of the voltage transfer functions. 
The decision was made to use both of these transfer functions. 

The factors leading to this decision are now given. 

In general, a two-port network would have one set of terminals 
for the normal input and the other set for the normal output. In 
this situation, the reverse voltage-transfer function may be con¬ 
sidered to be of very minor importance. However, it may provide 
data that is useful or even necessary to identify a faulty component. 

For any given network, the possibility exists that an element 
or elements of that network may not appear in the expression of 
one or more of the network functions. As an example, consider the 
ladder network of Figure 2. 
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Figure 2 


The voltage transfer functions, as defined by equations (1) 
and (2), can be found by a number of methods [3], [4], These have 
been calculated and are given in equations (3) and (4), 


R^R4 


’ Rl(R^-^R3+RJ +Rz(i?3+R-,) 


(3) 


R2 




(4) 


An inspection of these equations shows that elements and 
do not appear in the equation for For the proposed fault- 

identification procedure this means, of course, that no information 
can be obtained about the values of and R^ by making tests on 


Another reason for selecting both of the voltage transfer func¬ 
tions is the possibility that a change in the value of a network ele¬ 
ment may be offset by a change in one or more other elements, 
with the result that the magnitude of the transfer function will re¬ 
main unchanged. 
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This can be illustrated by again considering the network of 
Figure 2. Assuming that all of the element values are equal to 1, 
equations (5) and (6) are obtained by substituting these values into 
(3) and (4). 


G 


12 " 


G 


21 ■ 


1 ^ 

5 


1 ^ 

2 


(5) 

( 6 ) 


Now assume that R 2 has increased in value to 3, that has 
increased in value to 2, and that Rj^ and R^ have remained unchanged. 
The voltage transfer functions can again be evaluated by substituting 
these new network values into equations (3) and (4). The results 
of this are given in equations (7) and (8). 



I 

5 

3 

5 


(7) 

( 8 ) 


'A comparison of equations (5) and (7) shows tlmt the magnitude 
of has remained unchanged, although one of the network ele¬ 

ments has been increased in value by 300% and another by 200%. 
However, when equations (6) and (8) are compared, the change in 
magnitude of G 2 j indicates that a change has occurred in at least 
one of the network elements R 2 or R^. 

It is noted that no useful information for fault identification 
purposes would be obtained by making measurements of the voltage 
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transfer functions in the example given above, and it was stated 
earlier that purely resistive networks would not be allowed. However, 
this argument can readily be extended to cover the case of a network 
for which the magnitudes of its transfer fvmctions are frequency 
dependent. 

One other factor will be noted that argues for the use of both 
of the voltage transfer functions. It is, of course, the general 
catch-all phrase of having non-ideal test equipment and operators. 

In general, the test equipment should be operating within its per¬ 
formance specifications and the operator experienced enough in 
the use of this equipment so that the probability of an error in 
measurement would be lessened due to the larger number of measure¬ 
ments required from the use of both transfer functions. 

A measure of flexibility is added to the procedure through the 
use of both the forward and reverse voltage-transfer functions. 

The data for fault identification purposes will be in the form of two 
tables, one for each transfer function. It is possible that the in¬ 
formation needed can be determined by making measurements on 
only one of the transfer functions and the corresponding table for 
that transfer function. If the needed information cannot be deter¬ 
mined from this one set, the other would be available. It could 
also be used to serve as a check on the first set of measurements 
made. A form of this flexibility can be seen in the following 
example. 
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for this investigation. Therefore, equation (9) can be factored 
and given in the form of equation (10). 


H(s) 


(s- 


(10) 


The constant K is a gain term that is determined by the coefficients 
and b^^ of equation (9). The z and p terms are defined as the 
zeros and poles of the function. Although the coefficients of equation 
(9) are real, the poles and zeros of the function may be imaginary 
or complex. 

When imaginary or complex roots do occur, it is always in 
the form of conjugate pairs. This means that any imaginary or 
complex terms appearing in (10) may be eliminated by multiplying 
together the appropriate roots. When this is done, the network 
function is given in the form of equation (11). 


S*'(s- Pj)(s-Pt)- • • (s-Pj) • ' • ^ s V 6 ) 


( 11 ) 


Equation (11) is the form that is most useful for this study. 
The coefficients of all of the factors appearing in this equation are 
real. 


B. CORNER FREQUENCIES 

In this study the interest is in the sinusoidal steady-state 
operation of a network (i. e. , s = j ). Therefore, a frequency- 
response approach is used to investigate the factors making up 


18 











equation (11). Only the first-order and second-order factors are 
discussed, since they are the primary test frequencies. 

1. First-Order Factors 

Consider the first-order factor (s-z*) in equation (11). 

For sinusoidal steady-state operation, this factor can be re-defined 
by equation (12)* 


S- Z* £ JCO + We (12) 

As an example, consider the voltage transfer function given in 
equation (13). In order to use the frequency-response approach, 
the magnitude of this function can be expressed in dB and given by 
equation (14). 

Gj2 (j CO) = j£0 + 0)^ (13) 

At extremely low frequencies (60« (14) can be approxi¬ 

mated by (15). 

3iz~ ZOlo3„ll|:0 (15) 

For the extremely high frequencies ( 60 approxima¬ 

tion in (16) can be obtained. 

20lo3^Jj^| 
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The straight-line asymptotic approximation to (14) is given by 
(15) and (16). (16) is for the frequencies when and (16) is 

for frequencies where tO . A plot of this approximation is 

given in Figure 3, where the frequency axis is normalized with 
respect to C>>c. • 



The actual value of the transfer function at the frequency where 
CO = 60c, is 3 dB above the value of the straight-line approximation 
of Figure 3 (it would be 3 dBbelow the value of the straight-line 
approximation if the network function had consisted of a pole). The 
slope of the approximation is 0 for (0< 60^, and 20 dB per decade of 
frequency for CJ > 20 dB per decade of frequency when the factor 

is a pole). 

The frequency CO^ is defined as the corner frequency of a first- 
order factor. The corner frequencies for first-order factors are 
selected as some of the test frequencies. 
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2^ Second-Order Factors 

Tke second-order or quadratic factors of equatiotn (11) 
are considered. Factors of this form can be re-defined as in 
equation (17). 

= (Jaj)* + ? ]“60.(jc^)4 (On* (17> 

6 )n is defined as the undamped natural frequency, and ^ is 
the damping factor,. 

The quadratic factor can be investigated in the same manner 
that was used for the first-order factor. The transfer function 
considered, for this factor is given in equation (18),. 

^azCj^) ^ (j ( 18 ) 

Without carrying oxrt the steps for the investigation of the 
first-order factor, the results of the analysis on the quadratic 
factor are given in Figure 4 [T], 
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The slope of the curve in Figure 4 for CO < is 0; the slope 
is 40 dB per decade of frequency for 00>COr^ . The frequency axis 
has been normalized with respect to 60vi . 

The frequency 60^ is defined as the corner frequency for a 
quadratic factor. Frequencies resulting from this type of factors 
are used as test frequencies. 

Some comments on quadratic factors are now given. The 
quadratic factor has a peak value when the damping ratio is less 
than(VF^ )-l. The expressions for the magnitude of the peak 
value and the frequency at which it occurs are known [7], These 
are dependent upon the value of the damping ratio. Although these 
are important quantities and the peak value of the function can be 
large compared with the value of the straight-line approximation at 
that frequency, these have not been considered in this study. When 
this condition occurs in a complex system, a simple analytic ex¬ 
pression relating the corner frequency and the frequency at which 
the peak value occurs may not be available, 

C. ADDITIONAL TEST FREQUENCIES 

Although the poles and zeros of a function uniquely determine 
that function (to within an all-pass function) [l], the corner frequen¬ 
cies are not adequate in themselves to determine the necessary 
data for fault identification. Therefore, it became necessary to 
select additional test frequencies. 
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The selection of these additional test frequencies was done in 


a simple and arbitrary manner. However, very encouraging results 
have been obtained from the sample problems that have been worked 
using this method. 

The selection of these additional test frequencies was done in 
subroutine RESULT of the computer program, and this method is 
now described. 

The corner frequencies were calculated and placed into array 
WA. The array was then sorted nume rically by eliminating zero 
frequencies and placing the first non-zero frequency into the first 
position in the array. The zero frequencies were eliminated because 
of the sinusoidal steady-state operation being considered. The 
non-zero frequencies were compared, and when identical frequencies 
occurred, all but one were eliminated. 

Two additional frequencies were selected as follows: the lowest 
non-zero corner frequency was divided by ten, and the highest cor^ 
ner frequency was multiplied by two. 

The other frequencies were selected by comparing the ratio of 
adjacent corner frequencies. When this ratio was less than two, 
no additional frequency was chosen. When the ratio was less than 
ten, one additional frequency was selected by taking the arithmetic 
mean of the corner frequencies. When the ratio was less than 100, 
two additional frequencies were chosen as follows: one was twice 
the value of the lower corner frequency and the other was half the 
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value of the higher corner frequency. When the ratio was greater 
than 100, three additional frequencies were selected as follows: 
one was ten times the value of the lower corner frequency, one was 
one-tenth the value of the higher corner frequency, and the third 
was the arithmetic mean of the corner frequencies. 

All of these additional frequencies plus the corner frequencies 
were placed into array WB. This array was then sorted in the same 
manner as WA, The ratios of the adjacent frequencies in array 
WB were compared, and the higher of the two frequencies were elim¬ 
inated if this ratio was found to be less than 1. 01. The frequencies 
then remaining in array WB constituted the set of test frequencies. 

D. OTHER POSSIBLE TEST FREQUENCIES 

The selection of the additional test frequencies, as given in 
the previous section, was completely arbitrary. Although the 
results obtained using this method have been very good, another 
completely arbitrary method co\ild have been used. Perhaps fewer 
test frequencies would be required, and the results obtained might 
be as good or better. 

Another possible source of test frequencies would be to inves- 

I 

tigate the frequency response of the function to obtain the frequencies 
at which relative extrema occur [8]. At frequencies near these 
extrema, the slope of the curve should not be changing very rapidly. 

This would mean that an error in the setting of the frequency of a 
signal generator would not be as probable to give misleading information. 
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IV. FAULT ISOLATION SIGNATURES 


A. GENERAL 

Once the test frequencies have been selected, the next step 
is to use these frequencies to obtain the necessary data for fault 
identification purposes. 

The magnitudes of the voltage transfer functions are calculated 
for the nominal network values at each of the test frequencies. 

These magnitude ratios are then expressed in dB by multiplying 20 
times the common logarithm of the magnitude ratio. 

Next, the network elements have to be varied in value. This 
can be done over whatever range and in any combinations that are 
desired. The transfer functions are calculated for each variation 
and evaluated at the same set of test frequencies (also expressed 
in dB). 

The data obtained from these steps can be handled in any way 
that provides desired information. 

B. A SET OF SIGNATURES 

The evaluations of the transfer functions expressed in dB (or 
as a simple ratio) as given above would not be extremely useful in 
this form. For example, if this data were put into the form of a 
table to be used for comparison purposes, the table would be lengthy 
and difficult to read. Also, it is doubtful that an exact match would 
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be obtained when measurements were made on an actual network 
with the entries of the table. 

Therefore, it was decided to associate a single integer number 
(signature) for each value of the transfer functions at each of the 
test frequencies. This results in an n by m table where m is the number 
of test frequencies and n is the number of variations made on the 
network elements, and each entry in the table is a single integer 
number. 

The signatures were determined in the following manner. 

First, the gains of the functions for the varied elements were com¬ 
pared to the gains of the nominal network at the test frequencies. 

This was done by taking the absolute value of the difference between 
the nominal and varied gains. Next, signatures were assigned 
according to the format of Table I. 


Absolute value 
of difference 
in gain (dB) 


Signature 

assigned 


less than 0. 5 


0 

1 

2 

3 

4 

5 

6 

7 

8 
9 


0. 5 to 1. 0 
1.0 to 1. 5 

1. 5 to 2. 0 

2. 0 to 3. 0 
3; 0 to 4. 0 
4. 0 to 5. 0 
5.0 to 6. 0 
6. 0 to 8. 0 


greater than 8. 0 


Table I 
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The assigning of these signatures has been arbitrary. H only 
a specific type of network were to be studied, a specific set of 
signatures could be devised for this particular problem. At least 
one different, but similar method to the one given here has been 
used [l]. In this method gains above and below the nominal value 
were considered separately, and the signatures were assigned 
accordingly. 

An ’’arbitrary" method of assigning signatures may be 
necessary, because different types of networks may be assigned 
signatures from the same program. However, for a complex 
system it may be necessary to develop a specific set of signatures^ 
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V. RESULTS OBTAINED 


A. DISCUSSION OF THE PROGRAM 

A brief discussion of the program that is presented is now 
given. As given in Chapter II, the program taken from [3] was 
modified to generate only the voltage transfer functions. The only 
major changes made to this program were in MAIN and subroutine 
RESULT. Therefore, a discussion of the operation of the remaining 
subroutines is not given. 

Essentially the operation of the program up to the point where 
RESULT is called was not modified. When RESULT is called, the 
transfer functions for the nominal network are calculated in poly¬ 
nomial form. These are printed as output data. Then, the IBM 
supplied subroutine POLRT is called to evaluate the poles and zeros 
of the transfer functions. These roots are used to determine the 
set of test frequencies which are stored in array WC. The transfer 
functions for the nominal network are then evaluated at the test 
frequencies. The test frequencies (in Hz) and the corresponding 
values of the transfer functions (in dB) are printed as output data. 

The computation is then returned to MAIN. 

In MAIN a counter, LA, is updated and one of the network 
elements is varied in value. The tree-finding portions of the program 
are again repeated until RESULT is called. In RESULT the transfer, 
functions for the modified network are calculated in polynomial form. 


28 





These are evaluated at the test frequencies and s^ignatures assigned 
by comparing these values with those far the nominal network- The 
following is then printed as output: the element varied^ its deviiatian 
from the nominal value, and the associated signatures- This pro?- 
cess is repeated for each variation of an element- 

The elements were varied near the. end of MAIN, and only 
one element was varied at a time. By placing nested DO locrps after 
statement number 210 in MAIN, the elements can be varied in wfeaEt- 
ever combination is desired.. However, this was not done due to 
the excessive amount of printed output and computing time that 
would have been required, 

B-. EXAMPLES 

Two examples have been included that are typical of the results 
that have been obtained- These results are given in the COMPUTER 
OUTPUT section. 

1 • Example 1 

The network used for Example 1 is shown in Figure 5- 



Z9 
























In Figure 5 the nominal element values have been labelled. 

The element number has been given in a circle by its corresponding 
element, and the node numbers are given as required in [3], 

Z. Example 2 

The network used for Example 2 is given in Figure 6. 

The same discussion given in Example 1 applies here. 



The results obtained from these examples and others have 
been very encouraging. It can be seen from the COMPUTER OUTPUT 

section that the varied element can be easily identified. 
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VI. CONCLUSIONS 


The results obtained from this study indicate that the method 
presented here should be of practical value in a maintenance pro¬ 
gram (either preventive or corrective)- The method of selecting 
test frequencies has been shown to be a sound one- The assignment 
of signatures can be made flexible to meet various requirements by 
having a few interchangeable sets of signatures to be entered into 
the computer program. 

The computer program is not practical in its present form- 
It should be changed to generate the transfer functions in symbolic 
form. This is necessary to insure a reasonable computing time 
when the large number of combinations for the variations are con¬ 
sidered, because the tree-finding portions of the program could 
then be bypassed. Also, modifications need to be made to allow 
for the inclusion of mutual inductances and active network elements. 

The output of the program was left in a simple form because 
of the variety of ways in which it can be used. For example, the 
set of signatures (with the associated elements) could be sorted 
into some logical sequence to be placed into a table that would be 
used to look up faulty elements from actual measurements made on 
a circuit. Also, by expanding the program to allow the measurements 
made on a circuit to be used as inputs, the computer coiild compare 
these with its calculated set of signatures to provide only the faulty 


elements as output. 
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NODES= 5 


EDGE NUMBER 
1 
2 
3 
A 
5 

o, 6 

ro 

7 

A MATRIX 
10 10 
0 10 1 
0 0 0 0 
1111 
0 0 0 0 


i 


INPUT DATA 


EDGES= 7 GROUND CONDITICN= 


NODE NUMBERS 

1 4 

2 4 

1 4 

2 4 

4 5 

3 5 

2 3 


ADMITTANCE 

7.0000 

9.0000 

1.0000 

2.0000 

3.0000 

5.0000 

4.0000 


POWER 

1 

1 

0 

0 

0 

1 

0 


0 C 0 
0 0 1 
oil 
10 0 
lie 


i 




VOLTAGE TRANSFER FUNCTIONS 


6.00S** 0+ 79,00S** l+304,00S** 2+315.OOS** 3 
42.00S^* 0+421.OOS** 1+759.OOS** 2+315.OOS** 3 


TEST FREQUENCY IHZI 

GAIN fOB) 

0.203E-02 

-0.169E 02 

0.203E-01 

-0.161E 02 

0.227E-01 

-0.159E 02 

0.354E-01 

-0.146E 02 

0.654E-01 

-0.H4E 02 

C.955E-01 

-0.896E 01 

O.IOIE 00 

-0.858E 01 

0.182E 00 

-0.483E 01 

0.262E 00 

-0.299E 01 

0.524E 00 

-0.964E 00 


6.00S** 0+ 79.00S** 1+304.OOS** 2+315.OOS** 3 
6.00S** 0+ 94.00S** 1+409.OOS** 2+315.OOS** 3 


TEST FREQUENCY IHZ) 

GAIN (OB) 

0.207E-02 

-0.270E-01 

C.207E-01 

-0.160E 01 

0.227E-01 

-0.178E 01 

0.354E-01 

-0.257E 01 

0.654E-01 

-0.293E 01 

0.955E-01 

-0.257E Cl 

0.163E 00 

-0.160E 01 

0.326E 00 

-0.579E 00 








FAULT ISOLATICN SIGNATURES FCR V2/V1 


ELEMENT DEVIATION 


SIGNATURES 


CO 


1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

3 

3 

3 

3 

‘3 

3 

3 

3 

3 

3 

3 

3 

3 


-0.9 

-0.7 

-0.5 

-0.3 

- 0.1 

0.1 

0.3 

0.5 

0.7 

0.9 

1.1 

1.3 

1.5 

-0.9 

-0.7 

-0.5 

-0.3 

- 0.1 

0.1 

0.3 

0.5 

0.7 

0.9 

1.1 

1.3 

1.5 

-0.9 

-0.7 

-0.5 

-0.3 

- 0.1 

0.1 

0.3 

0.5 

0.7 

0.9 

1.1 

1.3 

1.5 


0 2 
0 2 
0 2 
0 1 
0 0 
0 0 
0 1 
0 2 
0 3 

0 3 

0 4 

0 4 

0 4 

0 0 
0 0 
0 0 
0 0 
0 0 
0 0 
0 0 
0 0 


0 1 
0 2 
0 2 
9 3 

8 3 

5 2 

3 1 

1 0 
0 0 
2 1 

3 2 

4 3 

4 4 

5 4 

5 4 


3 5 8 9 

3 4 7 8 

2 4 5 5 

12 3 3 

0 0 11 

0 0 0 0 

12 2 2 

2 3 3 3 

3 4 4 4 

4 4 4 4 

4 4 5 4 

4 5 5 5 

5 5 5 5 

0 2 5 6 

0 2 4 5 

0 13 4 

0 12 2 

0 0 0 0 

0 0 0 0 

0 111 

112 2 
12 2 2 

12 3 2 

2 3 3 3 

2 3 4 3 

2 4 4 3 

2 0 0 0 

2 0 0 0 

2 0 0 0 

10 0 0 

0 0 0 0 

0 0 0 0 

10 0 0 

2 10 0 

3 2 0 0 

3 2 0 0 

4 3 10 

4 3 10 

4 4 10 


9 9 9 9 

8 8 7 5 

5 5 4 2 

3 3 2 1 

10 0 0 

0 0 0 0 

2 110 

3 2 10 

4 3 2 0 

4 3 2 0 

4 4 2 1 

5 4 3 1 

5 4 3 1 

6 8 8 7 

5 5 4 2 

4 3 2 1 

2 110 

0 0 0 0 

0 0 0 0 

10 0 0 

110 0 
2 10 0 

2 110 

3 110 

3 110 

3 2 10 

0 0 0 0 

0 0 0 0 

0 0 0 0 

0 0 0 0 

0 0 0 0 

0 0 0 0 

0 0 0 0 

0 0 0 0 

0 0 0 0 

0 0 0 0 

0 0 0 0 

0 0 0 0 

0 0 0 0 






ELEMENT 


u> 

Ol 


O'<><><><> c^ ovjiunjf\;t»jfJLnu!vno!UTUiuivji^^-P^^ 


f 


DEVIATION 


SIGNATURES 


-0,9 

-0.7 

-0.5 

-0.3 

- 0.1 

0.1 

0.3 

0.5 

8:1 

1.1 

1:1 

-0.9 

-ol? 

0.1 

0.3 

C.5 

0.7 

:! 

{:i 

*-0*9 
-0.7 
-0.5 
-0*3 
-"9*1 
1*1 
• 3 
f»5 

8 :? 
if 

-Sis 

:8:J 

0*3 

0*5 

0*7 

III 

III 


9 

6 

4 

2 

0 

0 

1 

1 

2 

2 

2 

3 

3 

0 

0 

0 

0 

0 

0 

0 

0 


8 

0 

0 

0 

0 

0 


0 

8 

0 

0 

0 

0 

9 

8 

6 

4 
1 
1 
3 

5 

6 
6 

I 

8 


4 

4 

3 

1 

0 

0 

1 

2 

2 

3 

3 

3 

4 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
2 
1 
1 

8 

0 

0 


3 

3 

2 

1 

0 

0 

1 

2 

2 

3 

3 

3 

4 
1 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
2 
2 
1 
0 
0 
0 


I 

I 

6 

4 

1 

1 

3 

4 

5 


1 

1 

1 

0 

0 

0 

0 

1 

2 

2 

3 

3 

3 

2 

1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

3 

2 

2 

1 

0 

0 

1 

1 


0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1 

1 

1 

4 

2 

1 

8 

0 

0 

0 

0 

0 

0 

0 

0 

4 

3 

2 

1 

0 

0 

0 

1 


I 1 


0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

4 

3 
1 
0 
0 
0 
0 
0 
0 
1 
1 
1 
1 

4 
2 


0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

4 

3 

1 

1 

0 

0 

0 

0 

1 

1 


1 

4 

2 

1 

0 

I 

0 

0 

0 

0 

0 


? 


0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

3 


1 

0 

0 

0 

1 

1 

2 

2 

2 

2 

2 

1 

0 

0 

0 

0 

c 

0 


0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

2 

1 

1 

0 

0 

0 

0 

1 

1 

2 

2 

2 

3 

1 

0 

0 

0 

0 

0 

0 

0 


3 

4 

5 

I I 

8 8 




4 

5 
5 

I 


0 
0 
0 

4 
4 

I I 
I 8 

4 3 

I f 

1 6 
7 6 


8 8 


0 

0 

0 

3 

3 

2 

1 

0 

0 

I 

I 


0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1 


2 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

f 

0 

0 

0 

8 


3 

3 





FAULT ISOLATION SIGNATURES FOR V1/V2 
ELEMENT DEVIATION SIGNATURES 


O' 


1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 


-0.9 

-0.7 

-0.5 

-0.3 

- 0.1 

0.1 

0.3 

0.5 

0.7 

0.9 

1.1 

1.3 

1.5 

-0.9 

-0.7 

-0.5 

-0.3 

- 0.1 

0.1 

0.3 

0.5 

0.7 

0.9 

1.1 

1.3 

1.5 

-0.9 

-0.7 

-0.5 

-0.3 

-Q.l 

0.1 

0.3 

0.5 

0.7 

0.9 

1.1 

1.3 

1.5 


0 0 0 0 
0 0 0 0 
0 0 0 0 
0 0 0 0 
0 0 0 0 
0 0 0 0 
0 0 0 0 
0 0 0 0 
0 0 0 0 
0 0 0 0 
0 0 0 0 
0 0 0 0 
0 0 0 0 
0 110 
0 0 0 0 
0 0 0 0 
0 0 0 0 
0 0 0 0 
0 0 0 0 
0 0 0 0 
0 0 0 0 
0 0 0 1 
0 0 0 1 
0 0 0 1 
0 0 0 1 
0 0 0 2 
0 0 0 0 
0 0 0 0 
0 0 0 0 
0 0 0 0 
0 0 0 0 
0 0 0 0 
0 0 0 0 
0 0 0 0 
0 0 0 0 
0 0 0 0 
0 0 0 0 
0 0 0 0 
0 0 0 0 


0 0 
0 0 
0 0 
0 0 
0 0 
0 0 
0 0 
0 0 
0 0 
0 0 
0 0 
0 0 
0 0 
3 5 

3 4 

2 3 
1 1 
0 0 
0 0 
1 1 
1 1 
2 2 
2 2 
2 2 
2 2 

3 3 
0 0 
0 0 
0 0 
0 0 
0 0 
0 0 
0 0 
0 0 
0 0 
0 0 
0 0 
0 0 
0 0 


0 0 
0 0 
0 0 
0 0 
0 0 
0 0 
0 0 
0 0 
0 0 
0 0 
0 0 
0 0 
0 0 
7 7 

5 4 

3 1 

I 0 
0 0 
0 0 
0 0 


2 0 
0 0 
0 0 
0 0 
0 0 
0 0 
0 0 
0 0 
0 0 
0 0 
0 0 
0 0 
0 0 
0 0 






0 


ELEMENT OEVIATICN 


SIGNATURES 


-j 


4 

4 

4 

4 

4 

4 

4 

4 

A 

4 

4 

4 

4 

5 
5 
5 

c 

5 

5 

5 

5 

5 

5 

C 

5 

*; 

6 
6 
6 
6 
6 
6 
6 
6 
6 
6 
6 
6 
6 
7 
7 
7 
7 
7 
7 
7 
7 
7 
7 
7 
7 
7 


“0.9 

-0.7 

-0.5 

-0.3 

- 0.1 

0.1 

0.3 

0.5 

0.7 

0.9 

1.1 

1.3 

1.5 

-0.9 

-0.7 

-0.5 

-0.3 

- 0.1 

0.1 

0.3 

0.5 

0.7 

0.9 

I.l 

1.3 

1.5 

-0.9 

-0.7 

-0.5 

-0.3 

- 0.1 

0.1 

0.3 

0.5 

0.7 

0.9 

1.1 
1.3 
1.5 

-0.9 
-0.7 
-0.5 
-0.3 
- 0.1 
0 . 1 
0.3 
0.5 
0.7 
0.9 

1.1 
1.3 
1.5 


2 4 
0 3 
0 2 
0 1 
0 0 
0 0 
0 0 
0 1 
0 1 
0 1 
0 1 
0 2 
0 2 
0 1 
0 0 
0 0 
0 0 
0 0 
0 0 
0 0 
0 0 
0 0 
0 0 
0 0 
0 0 
0 0 
0 2 
0 2 
0 1 
0 1 


0 1 
0 2 
0 2 
0 3 

0 4 

0 4 

0 0 
0 0 
0 0 
0 0 
0 0 
0 0 
0 0 
0 0 
0 0 
0 0 
0 0 
0 0 
0 0 


3 

3 

2 

1 

0 

0 

0 

1 

1 

1 

2 

2 

2 

1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

3 

2 

1 

1 

0 

0 

1 

1 

2 

3 

3 

4 

4 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 


1 

1 

1 

0 

0 

0 

0 

1 

1 

2 

2 

2 

2 

3 
1 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 

4 
3 
2 

1 
0 
0 
1 
1 

2 
3 
3 

3 

4 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 


0 

0 

0 

0 

0 

0 

0 

0 

0 

1 

1 

1 

1 

4 

3 
1 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 

4 
3 
2 
1 
0 
0 
0 
1 
1 
1 
1 
1 
1 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 


0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

4 

3 
2 
1 
0 
0 
0 
1 
1 
1 
1 
1 
1 

4 
2 
1 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 


0 0 
0 0 
0 0 
0 0 
0 0 
0 0 
0 0 
0 0 
0 0 
0 0 
0 0 
0 0 
0 0 
3 1 

2 0 
1 0 
1 0 
0 0 
0 0 
0 0 
1 0 


2 1 
2 2 
2 2 
2 0 
0 0 
0 0 
0 0 
0 0 
0 0 
0 0 
0 0 
0 0 
0 0 
0 0 
0 0 
0 0 
0 0 
0 0 
0 0 
0 0 
0 0 
0 0 
0 0 
0 0 
0 0 
0 0 
0 0 
0 0 
0 0 


TIME= 30.31 SEC 






NODES* 4 


EDGE NUMBER 
1 
2 

3 

4 

5 

A MATRIX 
1 1 C 0 

10 10 
0 0 11 
0 10 1 




INPUT DATA 


EDGES* 5 GROUND CONDITION* 2 


NODE NUMBERS 
1 2 

1 4 

2 3 

3 4 

2 4 


ADMITTANCE 

1.0000 

1.0000 

2.0000 

1.0000 

1.0000 


PCWER 

1 

0 

-1 

0 

1 


0 

1 

0 

1 




VOLTAGE TRANSFER FUNCTIONS 


-2.00S** 0+ 

V2/V1 = - 

4.00S** 0+ 

TEST FREQUENCY (HZ) 

0,164E-01 
0.164E 00 
0.225E 00 
0.411E 00 
0,597E 00 
0,119E 01 

-2.00S** 0* 

V1/V2 = - 

2.00S** 0> 

TEST FREQUENCY (HZ) 

0,l59E-0l 
0.159E 00 
0.225E 00 
0.318E 00 
0.637E 00 


l.OOS** 2 

2.00S** 1+ 4ToOS**”2+ 

GAIN (DB) 

-0.589E 01 
0.973E 01 
0,249E-04 
-0.943E 01 
-0.127E 02 
-0.180E 02 

l.OOS** 2 

3.00S** T* I.oos** 2 ” 
GAIN (DB) 

-0.107E-01 
-0.458E 00 
-0.512E 00 
-0.458E 00 
-0.209E 00 


l.OOS** 3 







FAULT ISOLATION SIGNATURES FOR V2/VI 


ELEMENT DEVIATION 


SIGNATURES 


o 


1 

1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 

2 
2 
2 
2 
2 
2 
2 
2 
2 
2 
2 
2 
2 
3 
3 
3 
3 
3 
3 
3 
3 
3 
3 
3 
3 
3 


-0.9 

-0.7 

-0,5 

-0.3 

- 0,1 

0,1 

0.3 

0,5 

0.7 

0.9 

I.l 

1.3 

1.5 

-0.9 

-0.7 

-0.5 

-0.3 

- 0.1 

0.1 

0.3 

0.5 

0.7 

0.9 

1.1 
1.3 
1.5 

-0.9 

-0.7 

-0.5 

-0.3 

- 0.1 

0.1 

0.3 

0.5 

0.7 

0.9 

1.1 
1.3 
1.5 


0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

9 

8 

5 

3 

0 

0 

2 

3 

3 

4 
4 

4 

5 
2 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 


9 

9 

3 

5 

1 
0 
0 

2 

5 

6 

7 

8 
8 
8 
7 
5 

3 
1 
1 
2 

4 

5 

5 

6 
5 
7 
9 
9 
9 
7 
1 
1 
5 

7 

8 
9 
9 
9 
9 


2 

2 

3 

2 

0 

0 

2 

3 

4 

4 

5 
5 
5 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
9 
9 
8 
5 
2 
2 

5 
8 
8 
8 
7 

6 
5 


5 

2 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

2 

2 

1 

0 

0 

0 

0 

1 

1 

2 

2 

2 

3 
5 

4 
4 
2 
0 
0 
2 
4 

4 

5 

6 

7 

8 


7 

3 

1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

2 

1 

1 

0 

0 

0 

0 

1 

1 

2 

2 

3 

3 

3 

2 

2 

1 

0 

0 

1 

2 

2 

3 

4 

4 

5 


8 

3 

1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

.0 

0 

0 

0 

0 

0 

1 

1 

1 

2 

2 

1 

0 

0 

0 

0 

0 

0 

0 

0 

I 

1 

1 

1 





ELEMENT 


VJIOtUIVJiUlUlUIVJl V/IVTtUlOl 


$ 


DEVIATION 


SIGNATURES 


-0.9 

-0.7 

-0.5 

-0.3 

-O.l 

O.I 

0.3 

0.5 

0.7 

0.9 

1.1 

1.3 

1.5 

-0.9 

-0.7 

-0.5 

-0.3 

- 0.1 

0.1 

0.3 

0.5 

0.7 

0.9 

1.1 

1.3 

1.5 


7 

5 

4 

2 

0 

0 

2 

3 

4 

5 

5 

6 
6 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 


5 

4 

3 

1 

0 

0 

1 

2 

3 

4 
4 

4 

5 

6 
4 
1 

4 
2 
2 

5 

7 

8 
9 
9 
9 
9 


0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

9 

9 

8 

5 

1 

1 

4 

5 

6 

7 

8 
8 
8 


8 

6 

4 

2 

0 

0 

1 

2 

3 

3 

4 
4 
4 
8 
7 
«; 

3 

1 

1 

3 

4 

5 

6 
6 

7 

8 


9 

8 

5 

3 


2 

3 

4 

4 

5 
5 
5 
9 
7 
5 
4 
1 
1 

3 

4 

5 

6 
7 

7 

8 


9 

9 

7 

4 

1 

1 

3 

4 

5 

6 
7 

7 

8 
9 
8 
6 
4 
1 
1 

4 

5 

6 
7 

7 

8 
8 





FAULT ISOLATION SIGNATURES FOR VI/V2 


ELEMENT 


1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

2 

2 

2 

2 

2 

4^ 2 

N 2 

2 

2 

2 

2 

2 

2 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 


DEVIATION 


SIGNATURES 

-0.9 

0 

0 

1 

2 4 

-0.7 

0 

0 

0 

0 0 

-0.5 

0 

0 

1 

0 0 

-0.3 

0 

0 

0 

0 0 

-O.l 

0 

0 

0 

0 0 

0.1 

0 

0 

0 

0 0 

0.3 

0 

0 

0 

0 0 

0.5 

0 

1 

1 

1 0 

0.7 

0 

2 

1 

1 0 

0.9 

0 

2 

2 

1 0 

I.l 

0 

3 

2 

1 0 

1.3 

0 

3 

3 

2 0 

1.5 

0 

4 

3 

2 0 

-0.9 

4 

6 

5 

4 1 

-0.7 

0 

4 

4 

2 0 

-0.5 

0 

2 

2 

1 0 

-0.3 

0 

1 

1 

0 0 

-0.1 

0 

0 

0 

0 0 

0.1 

0 

0 

0 

0 0 

0.3 

0 

0 

0 

0 0 

0.5 

0 

0 

0 

0 0 

0.7 

0 

0 

0 

0 0 

0.9 

0 

0 

1 

0 0 

1.1 

0 

0 

1 

0 0 

1.3 

0 

0 

0 

0 0 

1.5 

0 

0 

0 

0 0 

-0.9 

1 

2 

1 

0 0 

-0.7 

0 

0 

0 

0 0 

-0.5 

0 

0 

1 

0 0 

-0.3 

-0.1 

0 

A 

0 

A 

0 

A 

0 0 

A A 

u 

0 

u 

0 

0 

0 

0 0 
0 0 

0.1 


n 

A 

0 

0 0 

U • <5 

U 

u 

0.5 

0 

1 

1 

1 0 

0.7 

0 

1 

1 

2 1 

0.9 

0 

1 

2 

2 1 

1.1 

0 

1 

2 

3 2 

1.3 

0 

2 

3 

3 3 

1.5 

0 

2 

3 

4 3 






ELEMENT 


4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

5 
5 
5 
5 


5 

5 

5 

5 

5 

5 

5 


TIME= 10.39 SEC 


DEVIATION 


SIGNATURES 
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0 0 0 0 
0 0 0 0 






INTEGER 

DIMENSION KAMI12,20),IS(20)tY(20»♦W(2C ),A(12»20),WV(9,20)»ISPI20)t 

1 ISN(20),LINE(100),0EN(20),onU(20),TITLE(2,9),NAME(6,10),YA(20) 

2 ,PCV(150),IE(150) 

COMMON ITF(2),ISIG(2,150,25),NTR=E(900,20) 

250 FORMAK/ 5X,6H TI ME=, P6.2 ,4H SEC) 

READ (5,49) BLANK,PLUS,DASH 
4<? FORMAT (3A1) 

DO 45 J=l,lOO 
45 LINE(J)=DASH 

READ (5,54) ((TITLE!I,J),J=l,8),1=1,2) 

54 FORMAT (16A1) 

READ (5,55) ((NAME(I,J),J=1,10),1=1,6) 

55 FORMAT (60A1) 

IICO CLCCK = ITIMF(0)4.01 

TIME=ITIME(0)♦.Ol-CLOCK 

IA=0 

C=-l.l 

JA=l 

READ (5,50) NI,MR,IGR,IPnRT 
50 F0RMAT(4T10) 

IF (MI.EQ.O) GO TO 1111 
PRINT 59 

59 FORMAT (IHl ,30X,'INPUT DATA'////) 

PRINT 7,NI,M8,IGR 
ITERM=2*NI>1 

4^ 7 FORMAT (/lOX,'NOOES = ',I2,5X,'EDG=S=', I2,5X,'GPOUND=', 12///, lOX, 

4^ I'EDGF NUMBER ' ,5X,'NODE NUMBERS ' ,5X, ' ADMITTANCE',5X,'POWER ' ^ ) 

DO 60 T=l,NI 
DO 60 J=1,MB 

60 A(I,J)=0 
JS=1 


DO 61 J= 
READ (5i 
51 FORMAT ( 
PRINT 6- 

= 1,MB 

,51) LE,NX,NY,Y(LE),IS(LE) 
(110,213,FIO.4,15) 
,LE,NX,NY,Y(LE),IS(LE) 

6 FORMAT ( 

[6X,I10,10X,2I6,2X,F12.4,6X, 

A(NX,LE1 

=1 

A(NY,LE] 

1=1 

YA( J)=Y 

r j) 

IF (ISd 

_E)) 62,61,61 


62 JS=JS+1 
61 CONTINUE 

100 FORMAT (//10X,'A MATRIX'/) 
WRITE (6,100) 

DO 101 1=1,NI 

101 WRITE(6,1C2) (A(I,J),J=1,MB) 

102 FORMAT (13X,25I3) 








Nl = l 
N2=2 

IF (IGR.FO.l) GO TO 205 

N3=3 

NA=A 

GO TO 210 
2C5 N3=3 
NA=3 
210 ?«F=6 

230 DO 20 L=ltMF 

GO TO (12,13,1A,15»16,17),1 

12 Tl=Nl 
I 2~N3 

CALL MODNET (11,12,A,NI,MB,KAM,N) 
GO TO 200 

13 Il=N2 
T 2*NA 

CALL MOONET (I1,I2,A,NI,MB,KAM,N) 
GO TO 200 
lA n=Nl 

CALL^MCDNET(I1,T2,A,NI,MB,KAM,N) 

GO TO 200 
I? I1=N3 
I2=NA 

IF (I1.E0.I2) GO TO A2 
CALL MODNFT (11,12 ,A, NT , MB, KAM, N) 
GO TO 200 
A2 N=1 

DO 72 I*1,ITERM 
72 WV(L,I)=0.0 
GO TO 20 

16 Il=Nl 
I2=N^ 

CALL MOONETdl ,I2,A,NI ,MB,KAM,N) 
GO TO 200 

17 Il=N2 
T2=N3 

CALL M00NFT(I1,I2,A,NI,MB,KAM,N) 
GO TO 200 
2CC CONTINUE 

CALL TREE (N,MB,KAM, IS , JS ,’V 

00 66 T=1,ITERM 

66 WV(L,I)=W(I) 

2C CONTINUE 

IF (lA) 68,69,68 
65 PRINT 67 

67 FORMAT (IHl,20X,'VOLTAGE TRANSFER 


,NI,W) 


FUNCTIONS'///) 






te on 86 I=1,ITKRM 

PNUd )=0.5*(WV(5,T»+WV(6,I )-WV(3,I)-WV(4, I) ) 

86 OENCI )=WV(l ,I ) 

LK=l 

CALL RESULT(DEN,PNU,NT,LINE,TSNf I?OtLK,TITLE, RL ANK , PLUS f I't JA ,C ) 
DO 87 I=1.TTERM 

PNUd )=C.5*(WV(5 ,I)+WV(6,n-WV(3,n-WV(4, I) ) 

87 DEN( I )=WV(2 ,n 
LK=2 

CALL RESULT(DEN,PNU,Ml,LINE,ISNtISD,LK,TTTLE,RLANK,PLUS,lA,JA,C) 
IA=IA+1 
GO rr 300 
3C2 C=-l.l 

300 C=C-l-0.2 

IF (C.GT.l.S) GO TO 301 
Y( JA)=YA( JA»+C*YA( JA) 

PCVdA)=C 
IEdA)=JA 
GO TO 230 

301 Y(JA)=YA(JA) 

JA=JA + 1 

IF (JA.GT.MB) GO TO 700 
GO TO 302 
700 IA=IA-1 

DO 701 J=l,2 

write (6,702) (TITLE!J,L),L=1,5) 

702 FORMAT (IHl,10X,•FAULT ISOLATION SIGNATURES FOP •,5A1///) 

WRITE (6,705) 

705 PORMAT dOX ,’ELEMENT* ,5X ,» DEV I AT ION*, lOX, • S IGNATURES*//) 

M=ITF(J) 

DO 703 L=1,IA 

703 WRITE (6,704) IE ( L ) , PCV ( L ) , ( IS IG ( J , L, K ), K= 1, M ) 

704 FORMAT (13X , 12 ,10X,F5.1,lOX,2513) 

TCI CONTINUE 

loco time=itimE(0)*.01-CL0CK 
WRITF(6,250) TIME 
GO TO 1100 
1111 STOP 
END 




4 


4 


4 


« 




SUBROUTINE MODNET (11♦I 2,A,NI,MB,KAM,N» 
INTEGER A 

DIMENSION A(12*20) ,KAM(12,20) 

DO 4 I =l,NI 
DO 4 J =l,MB 
4 KAH(I,J}sA(I,J) 

DO 3 J=1,MB 

KAMIIlfJ>=A(Il,J)+AII2,J) 

KAM(I2,J>=A(NI,JI 
IE (KAMdlfJ) .LE.ll GO TO 3 
KAM(I1,J)=0 
3 CONTINUE 
N=NI-1 
RETURN 
END 


it^ 


SUBROUTINE TREE (N.MB.KAM,IS,JS♦Y,NI,W) , ^ ^ ^ 

DIMENSION KEY(12) fKSa2,20 KKAMdZf 20),LTMP(12, 20),C0NTL( 12), 
1 IS(20),W(20),Y(20),KA(12,12,20) 

COMMON ITF(2),ISIG(2,150,25),NT REE(900,20) 

ITERM=24NI-1 
DO 58 K=l,ITERM 

58 W(K)=0.0 
NV=N-1 
KSUM=0 

00 59 1=1,NV 
DO 59 J=1,MB 

59 KA(1,I,J)=KAM(I+1,J) 

DO 60 J=1,MB 

60 KS(1,J)=KAM(1,J) 

NX1=N-1 

NX2=N-2 

NX3=N-3 

NX4=N-4 

NX5=N-5 

NX6=N-6 

NX7=N-7 

NX8=N-8 

NX9=N-9 

NX10=N~10 

NX11=N-11 

CONTL(1)=1.0 

DO 100 11=1,NXl 

NX=NX1 

M=1 






00 


IF (C0NTL( 
CALL 

IF (NX2.F0 
IF (C0NTL( 
DO 99 12= 
NX=NX2 
M«2 

IF (CONTLI 
CALL 

IF (NX3.E0 
IF (CONTLI 
DO 98 13= 
NX=NX3 
M=3 

IF (CONTL( 
CALL 

IF (NX4.E0 
IF (CONTLI 
DO 97 IA= 
NX=NX4 
M=4 


W),FO.O.O) GO TO 300 
NPART(M,M8,KEY,KS»KA,LTMP,C0NTL,NXI 
.0) GO TO 190 
M+1).EO.0.0) GO TO 300 
1 ,NX2 


M).EO.O.O> GO TO 100 
NPART(M,MB,KEr,KS,KA,LTMP,CONTL,NX) 
.0) GO TO 190 
M+n.EO.O.OJ GO TO 100 
1 ,NX3 


MI.EO.O.O) GO TO 99 
NOART(M,H8,KEY,KS,KA,LTMP,CONTL,NX) 
.0) GO TO 190 
M+H.EO.0.0) GO TO 99 
1 ,NX4 


IF (CONTL(M) .eO.0.0) GO TO 98 

CALL NPART(M,MB,KEY,KS,KA,LTMP,C0NTL,NX) 

IF (NX5.F0.0I GO TO 190 

IF (C0NTL(M+1).EO.O.O) GO TO 98 

DO 96 15 =1,NX5 

NX=NX5 

M=5 


IF (CONTL(M).EO.O.O) GO TO 97 

CALL NPART(M,MB,KEY,KS,KA,LTMP,CONTLtNX) 

IF (NX6.E0.0) GO TO 190 

IF (C0NTL(M+1).EO.O.O) GO TO 97 

DO 95 16 =1,NX6 

NX=NX6 

M=6 

IF (CONTL(M).EO.O.O) GO TO 96 

CALL NPARTCMfMB ,KEY t KS t KA ,LTMP, CONTL,NX ) 

IF (MXT.eO.O) GO TO 190 

IF (C0NTL(M+1).EO.O.O) 60 TO 96 

00 94 17 =1,NX7 

M=7 


IF (CONTL(W).EO.O.O) GO TO 95 
NX=NX7 

NPART<H,MR,KEY,KS,KA,LTMP,CONTL,NX) 
IF (NX8.E0.0) GO TO 190 
IF (C0NTL(M+1).EO.O.O) GO TO 95 
00 93 18 =1,NX8 
NX=NX8 


% 


i 


4 




4 


« 


vO 


190 

105 


205 

206 

87 

88 
89 

90 

91 

92 

93 
9A 


0 

IF (CONTLIWI.EO.O.O) go to 94 

CALL NOART(M,WB,KEY,KSfKA,LTMPtCONTL,NX| 

IF (NX9.EO.O) GO TO 190 
IF (CONTLCM+1).EO.O.O) GO TO 9^ 

00 92 19 =1,NX9 
NX=NX9 

IF (CONTLIM).EO.O.O) GO TO 93 , 

CALL NPARTCM,MBfKEY»KS,KA,LTMP,CONTL,NX) 

IFINXIO.EO.O) GO TO 190 

IF (C0NTL(M+1).EO.O.O) GO TO 93 

00 91 110=1tNXlO 

NX=NX10 

M=10 

IF (CONTLIM).EO.O.O) GC TO 92 

CALL NPART(M,HB,KEY,KStKA,LTMP,CONTL,NX) 

TFfNXll.EO.O) GO TO 190 

IF (C0NTL(M+1).EO.O.O) GO TO 92 

DO 90 111=1,NXll 

NX=NX11 

M=ll 

IF (CONTL(M).EO.O.O) GO TO 91 

CALL NPART(M,M6,KEY,KS,KA,LTMP,CONTL,NX) 

CONTINUE 

NTR=1 

DO 105 1=1,NV 
NTR=NTR*KEY(I) 

CALL XTREE(KEY,NTR,NV,MB,LTMP) 

DO 206 1=1,NTR 

V=1.0 

JJ=0 

DO 205 J=1,PB 

IF (NTREE<I,J).EO.O) GO TO 205 
V=V*Y(J) 

JJ=JJ*IS(J) 

CONTINUE 

W( JJ+JS)=V+WI JJ+JS) 

CONTI NUE 

GO TO (100,99,98,97,96,95,94,93,92,91,90,89,88,87),M 

CONTINUE 

CONTINUE 

CONTINUE 

CONTINUE 

CONTINUE 

CONTINUE 

CONTINUE 

CONTINUE 





95 CONTINUE 

96 continue 

97 CONTINUE 
95 CONTINUE 

CONTINUE 
ICO CONTINUE 
30C CONTINUE 
RETURN 
END 


SUBROUTINE NPART(M,M8,KEY,KS♦KA,LTMP,CONTL,NX) 
DIMENSION KEY(12),KS(12,20)#KA(12,12,20 I,LTMP(12,20), 
1 CCNTL(12) 

KEY( M)=0 

KEND=0 

KCPR=0 

NTEST=0 

DO 2 J=1,Me 

K=J 

IF (KS(M,J).EO.O) GO TO 2 
GO TO 9 

2 CONTINUE 

9 DO 3 1=1,NX 
IK=I 

TF(KA(M,I,K).EO.OIGO TO 3 
GO TO 15 

3 CONTINUE 

15 DO 72 J=1,MB 

LTMPf M,J)=KS(M,J)*KA(M,TK, J) 

KEY(M)=KEY(M)+LTMP(M, J) 

KA(M,IK,J)=KA(M,IK,JI-LTMP(M,J) 
KS(M,J)=KS(M,J)-LTMP(M,J) 

KS(M+1,J)=KS(M,J)+KA(M,IK,J) 

KEND=KEND+KS(M,J) 

KCPR=KCPR+KA(M.IK, J) 

NTEST=NTEST-i-KS{M+l ,J) 

72 CONTINUE 

DO 12 J=1,MB 
ITeMP = KA(M,NX, J) 

KA(M,NX,J)=KA(M,IK, J) 

12 KA(M,IK,J)=ITEMP 
K=NX-1 
DO 66 1=1,K 
DO 66 J=1,MB 
66 KA(M+1,I,J)=KA(M,I ,J) 




CQNTL(M)=KEND*KCPR 

CONTL(f!+l»=NTEST 

RETURN 

END 


20 


30 

28 

26 


SUBROUTINE XTREE < KEY,NTR,NV,MB t LTMP ) 
DIMENSION LTMP(12 ,201 .KEYdZ) 

COMMON ITF(2),ISIG(2,150,25),NTREE(ROO,20I 
DO 20 1*1,NTR 
DO 20 J=1,MB 
NTREEd ,JI=0 
MP = 1 

DO 26 1=1,NV 
MP=MP*KEY{I) 

MO=NTR/MP 
MS=MO*KEYd ) 

M=l 

DO 26 J=1,MB 

IF (LTMPd ,J)-1 >26,30,26 

DO 28 IM=M,NTR,MS 

DO 28 IJ=1,M0 

IS=IM4IJ-1 

NTREEdS,J)=l 

M=M+M0 

CONTINUE 

RETURN 

END 


SUBROUTINE RESULT(DEN,PNU,NI,LINE,ISN,ISO,LK,TITLE,BLANK,PLUS,lA, 

^OIMEnIiCN den( 20) ,PNU(20) ,LINFd0O>,ISN(2O), ISDt20),TITLE( 2,8), 

1 PSIGN(20),0SIGN(20),PPNU(20),PNUR(20),PNUI(20),ODEN(20), 

2 DENR(20),DENI(20),COF(20),WA(25),WB(2,25),B(20),H(2,25),HA(25), 

3 WC(2,25) 

COMMON ITF(2),ISIG(2,150,25),NTREE(900,20) 

I TERM=2*NI-l 
00 19 J=l,20 
PPNU( J)=0. 0 
19 ODEN(J)=0.0 
28 DO 29 J=1,ITERM 

IF (ABS(PNU(J)).GE.1000.) GO TO 30 
IF (ABS(DEN(J)).GE.1000.) GO TO 30 





29 CONTINUE 
GO TO 32 

30 DO 31 J=1,ITERM 
PNU(J)=PNU{J>/2. 
OEN(J)=DEN(J)/2. 

31 CONTINUE 
GO TO 28 

32 CONTINUE 
K=C 

I=C 

ESP=l.E-4 


DO 33 J =l,ITERM 

IF(ABS(DEN(JII.LT.ESP) GO TO 3 ^ 

K=K-H 

nEN(K)=DEN(J) 

ISD(K)=J 

IF (K.EO.l) GO TO 11 
IF (DENCK)) 11,12,12 

11 OSIGN(K|=BLANK 
GO TO 34 

12 09TGN(K)=PLUS 

34 IF(ABS(PNUU)).LT.ESP) GO TO 33 
1 = 1 + 1 

PNU(I»=PNU<J) 

ISN(T)=J 

IF (I.EQ.l) GC TO 18 
IF (FNUII)) 18tl7,17 
18 PS1GN( n=BLANK 
GO TO 33 

17 PSIGM(n=PLUS 
33 CONTINUE 

IF (ISDID-lSNd) ) 20,20t21 

20 M=isn(i) 

GO TO 25 

21 W=ISN(1) 

25 DO 36 J =1,K 

36 tSD( J)=ISDU)-M 

37 ISN?J)-M 

38 PRINT*72^n?llSN?J) t PNU (J) , IS N< JH t J= I I { 

PRINT 73t(TITLE!LK,J)♦Jflt9)♦(LINE(M),M=1, 
PRINT 74»(<QSIGN(J»,DEN{J),ISD(JJ)»J=ItK) 

72 F0RMAT(//20X,10(Al,F6.2,'S**'I2n 

73 FORMAT (10X,120AH 

74 FORMA TC20X,10(A1 ,F6 ,2 » ’ 12 )) 

n=isNn) 

KK=ISD{K) 


100 ) 




f 


a 


oi 


DO 301 J=1,I 
L=ISNCJ>+1 

301 PPNU(L)=PNU(J) 

00 302 J=ltK 
L=ISDtJ)+l 

302 DDEN(L>=0EN(J> , . 

CALL P0LRT(PPNU,C0F,II,PNUR,PNUI,lERN) 
CALL POLRT(DOEN,COF,KK,OENR,OENI,lERD) 
DO 501 J=1»II 

501 WA( J)=SQRT(PNUR( J» ♦♦2-»-PNUnj)**2l 
no 502 J=ltKK 

502 WA( J+II )=SORT(DENR( J)A>*2 + DENn J)**2) 
II*II+KK 

KK=n“l 
DO 601 J=1fII 
DO 602 P=1tKK 

IF CWA (W) . LE.WA( M-l-1) I GO TO 602 
ESP=WA(H+1) 

WA(H*1J=HA(M) 

WA(M)=ESP 
602 CONTINUE 
6C1 CONTINUE 

IF fWA(ll.EO.O.O) GO TO 803 
DO 801 J=1fKK 

IF (WA( J).NE.WAIJ4-1)) GO TO 801 
WAIJ)=0.0 

801 CONTINUE 
L=C 

DO 802 J*!,!! 

IF |WA(JI.EQ.O.O) GO TO 802 
L=L+1 

HB(LK,L)=WA(J) 

802 CONTINUE 
GO TO 806 

803 WB(LK,1)=0.0 
DO 804 J=ltKK 

IF (WA( J).NE.WA( J+D) GO TO 804 
WA(JJ=0.0 

804 CONTINUE 
L=1 

DO 805 J=lfII 

IF (WA(J).EQ.0.0) GO TO 805 
L=L+1 

WB(LK,L»=WA(J) 

805 CONTINUE 

806 CONTINUE 

w=0 

II=L-1 








00 101 J=1»I! 

IF (WB(LK,J).EO.O.OI GO TO 101 
IF (WBILKfJ+1)/W8(LK,J).LP.2.0» GOTO 101 
IF (WB(LK,J+1)/W8(LK»JI.LE.IO.C) GOTO 102 
IF (WB(LK,J+1)/W8(LK,J).LE.100.0) GOTH 1C3 
M=M+1 

(LK,L+M) =WB(LK , J) *10. 

M=Vi+l 

W8(LK,L+M) =WB(LK, J+1) /lO.O 
M=M+1 

WB(LK,L+f<) =(WB(LK,J) + WB(LK,J+1) )/2.0 
GO TO ICl 

103 M=M^1 

WB(LK,L+H)=WB(LK,J)*2.0 

M=M+1 

W8(LK,L+M)=WB(LK , J+1)/2.0 
GO TO 101 
102 M=M^1 

WB(LK,L+M)=(W8(LK,J)+WR(LK,J+1))/2,0 
ICl CONTINUE 
N=M+1 

WB(LK,L+M)=ViB(LKtl)/10.0 

M=M+1 

WB(LK,L+M)=WB(LK,L)*2.0 

L-L+M 

M=L-1 

DO 104 KK=1,L 
no 105 J=1,H 

IF ( WR(LK, J) .LE.WBILK ,J<-1) ) GOTO 105 
ESP=WB(LK,J+1) 

WB(LK,J+1)=WR(LK,J) 

WB(LK, J)=F5P 
1C5 CONTINUE 

104 CONTINUE 
M=C 

KK^L-l 

00 104 J=1fKK 

IF (WB(LK,J).EO.G.O) GO TO 106 

IF (WB(LK,J+1)/WB(LK,Jl.LT.l.Ol) GOTO 106 

M=M+1 

WCILKtM) =WB(LK,.J) 

106 CONTINUE 
M=M+1 

WC(LK,V)=WB(LK,L) 

ITF(LK)=M 
200 DO 210 J=l»20 
210 B(J)=0.0 
L=ITF(LK) 




« 




« 


A 

230 


Ol 

Oi 


8 

2A0 

221 

222 

220 


DO 220 M=1,L 
X«=0,0 
XI =0.0 
YR=0.0 
Y1=0.0 

B?JI=PNuTjf*WCaK,M)**ISN(J) 

IF nSNIJ).EO.O» GO TO A 

GO^TO^ (i f2f3fAflt2»3,4,lt2f3»4flt2f3»4*l»2»3)*N 

XI=XI+B(J) 

GO TO 230 
XR=XR-B( J) 

GO TO 230 
XI=XI-B( J) 

GO TO 230 
XR = XR+B( J» 

CONTINUE 

B?jf=DEN(j{*WC(LK,MI»AISD(Jl 
IF (ISD(JI.EO.O) GO TO 8 

GO^ToWs ,6,7,8,5 t6tT,8f5»6f7*8,5t6t7, 8,5,6t7)tN 

YI=Y1+6(J) 

GO TO 240 
YR=YR-B(J) 

GO TO 240 
YI=YI-B(J) 

GO TO 2A0 
YR*YR+B< J) 

CONTINUE 

HA(ll!=20?oiitLOG10(SORT(XRAA2+XI*A2)/SQRT(YRA*2+YI*A2)) 

H( LK%) =20. O^ALOGIO (SORT IXR*^2+X I**2) /SORT (YR**2+YI**2 ) » 

CONTINUE ^ 

IF (lA.EO.O) GO TO 270 
DO 40 J=lfN 

HA(J»=ABS(H(LK,J»-HA( JM 
IF (HA(J).LT.0.5) GO TO 41 
IF (HA(JI.LT.l.OI GO TO 42 
IF (HA{J>.LT.1.8) GO TO 43 
IF (HA(JI.LT.2.0j GO TO 44 
IF (HA(J).LT.3.0j GO TO 45 
IF (HA(JI.LT.A.O GO TO 46 
IF (HAIJ .LT.5.0) GO TO 47 
IF IHA(J).LT.6.0) GO TO 48 
IF (HA(J).LT.8.0) GO TO 49 







Ul 

o> 


4 P 


47 

46 

45 


44 

43 


42 

41 

40 

270 

271 

250 

260 

272 
2C'? 


I<;iG( LKtIA , J)=9 
GO TO 40 
ISIGf LK,IA,J)=8 
GO TO 40 
ISIG(LK,IA,J)=7 
GO TO 40 
ISIG(LK,IA,J)=6 
GO TO 40 
ISIGCLK,IA,J»=5 
GO TO 40 
ISIG(LK.IA,J)=4 
GO TO 40 
ISIG(LK,IA,J)=3 
GO TO 40 
ISIG(LK,IA,J)=2 
GO TO 40 
ISIG(LK.IA,J)=l 
GO TO 40 
ISIG(LK,lA tJ»=0 
CONTINUE 


GO TO 299 
00 271 J=1,M 

WC(LK,J)=WC(LK,J)/6.28318 
WRITE(6,250) 

F0RM4T(/10X.'TEST FREQUENCY {HZ)•,lOX,•GAIN 
PRINT 260,((WCfLK,J),H(LK,J)),J=1,M) 

FORMAT (14X,E10.3,15X,E10.3/) 

DO 272 J=1,M 

WC (LK,J)=WC(LK,J)*6.28318 

CONTINUE 

RETURN 

END 
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